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Abstract 

Using a lattice-gas description of the low-energy degrees of freedom of the quantum Heisenberg 
antiferromagnet on the frustrated two-leg ladder and bilayer lattices we examine the magnetization 
process at low temperatures for these spin models. In both cases the emergent discrete degrees of 
freedom implicate a close relation of the frustrated quantum Heisenberg antiferromagnet to the clas- 
sical lattice gas with finite nearest-neighbor repulsion or, equivalently, to the Ising antiferromagnet 
in a uniform magnetic field. Using this relation we obtain analytical results for thermodynamically 
large systems in the one-dimensional case. In the two-dimensional case we perform classical Monte 
Carlo simulations for systems of up to 100 x 100 sites. 

PACS numbers: 75.10.Jm 
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I. INTRODUCTORY REMARKS 



Experimental and theoretical studies of magnetization processes is a hot topic of the mod- 
ern condensed matter physics. Although a common wisdom says that the magnetization M 
monotonically increases approaching its saturation value as the applied field h increases, 
experimental observations for many materials often demonstrate a nontrivial dependence 
of M on h at low temperatures, which may show plateaus and jumps (see, e.g., Ref. 
This kind of behavior may be caused by competing interactions which occur in frustrated 
quantum systems. The theoretical description of the magnetization processes in frustrated 
quantum spin systems has attracted a lot of interest during the last decades, and theoret- 
ical predictions for the low-temperature magnetization curves of thermodynamically large 
systems are of great interest jlj]. 

On the other hand, recently the concept of independent localized-magnon states for the 
quantum Heisenberg antiferromagnet on certain classes of frustrated lattices has been in- 
troduced Q]. Later on it has been successfully used to describe magnetothermodynamics of 
these systems The independent localized-magnon states are the ground states in strong 
magnetic fields and therefore the independent localized-magnon picture (or hard-core model 
description) is adequate at the strong-field low-temperature regime. Moreover, very recently, 
the independent localized-magnon description has been improved: For some lattices (4] it 
appears possible to take into account the low-energy excited states (interacting localized- 
magnon states) too |5|. As a result, a lattice-gas description of low-energy degrees of freedom 
of certain frustrated quantum spin systems has been elaborated. In Ref. [5j some results 
for thermodynamic quantities (mainly concerning temperature dependences of specific heat, 
staggered susceptibility and entropy) have been reported too. The aim of the present paper 
is to use the lattice-gas description to study magnetization processes at low temperatures 
which have not been considered so far. We shall demonstrate that such a description for the 
considered frustrated quantum magnets yields very accurate results for the low-temperature 
quantities. 

The paper is organized as follows. In the next section (Sec. (TTJ) we define the quantum 
spin models and explain how the lattice-gas picture emerges. Then, in Sec. 1111} we calculate 
the magnetization of the frustrated quantum spin systems considering separately the one- 
dimensional case and the two-dimensional case. We report results for thermodynamically 
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large systems obtained using transfer-matrix method (one-dimensional case) and classical 
Monte Carlo simulations (two-dimensional case). We discuss in detail the low-temperature 
uniform magnetization and susceptibility within the frames of the lattice-gas description. 
Finally, in Sec. IIV[ we summarize our findings. 



II. FRUSTRATED QUANTUM SPIN SYSTEMS. LATTICE-GAS DESCRIPTION 

We consider N = 27V quantum spins 1/2 placed on the two lattices shown in Fig. [H see 



also Refs. 



7|. The Hamiltonian of the model reads 



H = ^J pq B p -s q -hS z , (2.1) 

(p?) 

where the sum runs over the bonds connecting neighboring sites on the lattice, J pq > is 
the antiferromagnetic interaction between the sites p and q, J pq takes two values, J2 for 
vertical bonds and J\ for horizontal and diagonal bonds, h > is the external magnetic 
field, S z = Ylp=i s p i s the z component of the total spin of the system. In our study we 
impose periodic boundary conditions and usually set J\ — 1 to fix the units. 

It is useful to introduce the total spin on a vertical bond t m = s m l + s m 2 . Here m = 
1 , . . . , M enumerates vertical bonds and runs over lattice sites of a simple chain in the case 
of the frustrated two-leg ladder or a square lattice in the case of the frustrated bilayer. Then 
the Hamiltonian (12. ip becomes 

H = E [y (*» - 1) - htz m] +*£w tl, (2.2) 

m L ^ ' -I (ml) 

where the second sum runs over the nearest-neighbor bonds on the simple chain or the 
square lattice. It is evident from Eq. ( 12. 2D . that the total spin of each vertical bond t m , 
m = 1, . . . , Af, t m = t m (t m + 1), is a good quantum number. Hence, the Hamiltonian (12. 2p 
depends on the set of quantum numbers {t m }, t m = 0, 1. As a result, the properties of the 
model can be studied in much more detail. In particular, a subset of 2^ = 1.41 . . . N low-lying 
eigenstates of the total set of 2 N eigenstates can be constructed exactly a a. Moreover, 
their contribution to thermodynamics can be calculated using a classical lattice-gas model 

fl. 

We start with an illustration of the eigenstates which we will consider in more detail. 
First we consider states which consist of n singlets on the vertical bonds m 1; . . . ,m n (i.e., 




FIG. 1: The two lattices considered in the paper: the frustrated two-leg ladder (a) and the frus- 
trated bilayer (b). We also show the phase diagrams as they follow from the lattice-gas description 
of the s = 1/2 Heisenberg antiferromagnet in a magnetic field on the frustrated two- leg ladder 
lattice with J\ = 1, Ji = 4 (h,2 = 4, hi = 6) and on the frustrated bilayer lattice with J\ = 1, 
J2 = 5 (/12 = 5, h\ = 9). A staggered occupation of vertical bonds by localized magnons (a kind 
of "antiferromagnetic" long-range order) occurs along the red line between /12 and hi at T = 
(ladder) or below the red critical line T c (h) with starting and end points at /12, T = and hi, T = 
(bilayer), whereas a uniform occupation of vertical bonds (a kind of "ferromagnetic" long-range 
order) occurs along the green lines h < /12 at T = and hi < h at T = 0. The remaining part of the 
phase diagrams corresponds to a disordered phase. Moreover, we show the lines below which the 
exact diagonalization data and the lattice-gas predictions for the specific heat of the finite system 
of iV = 16 sites coincide with the accuracy up to 5%, 2%, and 1% (blue dashed lines from top to 
bottom). 
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t^. = 0) and M—n fully polarized triplets (t^ = 2, t z m = 1) on the remaining vertical bonds. 
We impose a hard-core rule, i.e., the occupation of neighboring vertical bonds by singlets 
is forbidden. Following the notations introduced in Refs. js — 4 1 we call these states "inde- 
pendent localized-magnon states" and a singlet on a vertical bond a "localized magnon". 
The energy of these independent localized-magnon states is E l ™ = -Efm — n>£i, where 
Epm = N Ji+N ^2/4, ei = J 2 + 2Ji = hi (one-dimensional case) and Efm = 2A/ - J\+N "J2/4, 
ei = J2 + 4 Ji = hi (two-dimensional case). The degeneracy of the independent localized- 
magnon states g_\f(n) equals to the canonical partition function Zh c (n,N) of n hard-core 
objects (hard dimers or hard squares) on the lattice (simple chain or square lattice) of M 
sites. It is important to mention that the independent localized-magnon states are the 
ground states in the subspaces with S z = M — 1, . . . ,jV/2 if J2/J1 > 2 (one-dimensional 
case) or J2/ J\> 4 (two-dimensional case), cf. Ref. |8j, and that these states are linearly 
independent [9]. 



Next we relax the hard-core rule, i.e., we allow the localized magnons (singlets on vertical 
bonds) to be nearest neighbors. According to Ref. jsj] we call these states "interacting 
localized-magnon states". Any pair of neighboring localized magnons increases the energy 
by J\. As a result, for a given S z = M — n, n = 1, . . . , J\f/2 the energies of these interacting 
localized-magnon excited states become = M m + vJ\, where v is the number of pairs of 
neighboring localized magnons, cf. also Ref. (5j. The interacting localized-magnon states 
are the low-lying excited states for S z = Af — n, n = 1, . . . ,J\f/2 in the strong-coupling 
regime, i.e., when J2/J1 is sufficiently large. More precisely, from exact diagonalization data 
for finite systems we have found that the strong-coupling regime occurs when Jo > J 2 with 
J2/J1 ~ 3.00 (one-dimensional case) and J2/J1 ~ 4.65 (two-dimensional case) [5|. 

For lower values of the magnetization, S z = J\f/2 — r, r = 1, . . . ,J\f/2, where no indepen- 
dent localized-magnon states exist, the class of interacting localized-magnon states contains 
the ground-state manifold as well as low-lying excited states in the strong-coupling regime. 
The ground-state manifold is built by n = J\f/2 + r, r = 1, . . . ,J\f/2 localized magnons, 
where, e.g., Af/2 magnons occupy one sublattice of the underlying lattice (simple chain 
or the square lattice) completely, and the remaining r localized magnons sit on the other 
sublattice. The energy of this state is Ej^/ 2 + r = —MJ^I^ — rJ 2 , J2 = h 2 . The low-lying 
excited states are constructed from the ground state by rearranging the localized magnons to 
increase the number of neighboring magnons. Then each new pair of neighboring localized 



5 



magnons increases the energy by J\. 

Note that the interacting localized-magnon states can be visualized as partially overlap- 
ping hard-core objects (in contrast to the independent localized-magnon states which can be 
visualized as nonoverlapping hard-core objects). Although the degeneracy of the interacting 
localized-magnon states can be also calculated in terms of a canonical partition function of 
a system of hard-core objects, the required contribution of the independent and interacting 
localized-magnon states to the partition function Z(T, h, N) of the quantum spin system 
(12. ip is conveniently taken into account within the frames of a lattice-gas model of classical 
particles with finite nearest-neighbor repulsion V — Ji, see Ref. ij]. 

In the strong-coupling regime, when the constructed independent and interacting 



localized-magnon states dominate the partition function Z(TJi, N) of the quantum spin 
system (12.11) at low temperatures, for Z(T, h, N) we can write [5] 

Z(T,h,N)^Z LM (T,h,N)= ••• E e T 

ni=0,l 71^=0,1 

= e-^^S lg (T, / u,AT), (2.3) 

where Si g (T, fi, J\f) = J2 ni =o l • • ■ Snv=o i e~ n ^ nm ^^ T is the grand-canonical partition func- 
tion of the lattice gas with the Hamiltonian "H({n m }) = —^J2m n m + ^iE(mi) n "' I *i, 
li — h\ — h. Introducing the on-site spin variables a m = ±1 according to the relations 
a m = 2n m — 1 and n m = (1 + cr m )/2 one arrives at the Hamiltonian of the antiferromagnetic 
Ising model in a uniform magnetic field. 

In Ref. I^] we have performed extensive studies of the relation between the frustrated 
quantum spin model (12.11) and the classical lattice-gas model with the Hamiltonian H({n m }) 
for finite systems (up to N = 32) to clarify to what extent the latter model can reproduce the 
thermodynamic properties of the former one in the strong-coupling regime. In particular, 
we have checked the energies and degeneracies of low- lying states for S z = N/2, . . . , 0, have 
compared the field and temperature dependences of various thermodynamic quantities, have 
estimated the temperature range until which the lattice-gas predictions are in an excellent 
agreement with exact diagonalization data. Thus, for the frustrated two-leg ladder with 
Ji = 1, J2 = 4 (characteristic fields are — 4 and h\ = 6) and for the frustrated bilayer 
with J\ = 1, J2 = 5 (characteristic fields are /i2 = 5 and hi = 9) we have found that the 
lattice-gas picture works perfectly well at least up to T = 0.5 (and even up to T = 1 if 
J 2 has larger value J 2 = 10). To supplement this conclusion of Ref. js| we perform a new 
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comparison of the exact diagonalization data and the lattice-gas predictions for the specific 
heat of the finite system of N = 16 sites. Our results are shown in the phase diagrams in 
Fig. [TJ We find that below the dashed lines in the phase diagrams in Fig. [T] both results 
coincide with the accuracy up to 5%, 2%, and 1% (from top to bottom). 

The considered frustrated quantum spin systems exhibit an interesting low-temperature 
behavior if the magnetic field h is between h,2 and h\ which is related to an ordering of 
localized magnons on the simple chain or the square lattice. Both these lattices are bipar- 
tite ones, i.e., they consist of two sublattices A and B and any neighboring sites always 
belong to the sublattice A and the sublattice B. In the two-dimensional case the filling of 
the sublattices by localized magnons realizes as an order-disorder phase transition, which 
has a direct analogy to the phase transition for the square-lattice Ising antiferromagnet in 
a uniform magnetic field and belongs to the two-dimensional Ising model universality class. 
There are many studies on the square-lattice Ising antiferromagnet in a uniform magnetic 
field, see, e.g., Refs. 10Ml3|. and we may borrow the existing knowledge to examine the 
two-dimensional frustrated bilayer quantum Heisenberg antiferromagnet. In particular, a 
phase diagram in the half-plane "magnetic field - temperature" of the square-lattice Ising 
antiferromagnet in a uniform magnetic field has been discussed and a critical line T c (h) which 
separates the antiferromagnetically ordered phase and the disordered phase has been calcu- 
lated using different approximations. This phase diagram in the context of the frustrated 
bilayer has been reproduced in our classical Monte Carlo simulations j^j, see Fig. [Ib. While 
crossing the critical line, the thermodynamic quantities for the frustrated bilayer exhibit 
singularities: The specific heat shows a logarithmic singularity, the staggered magnetization 
(which may play a role of the order parameter) decays within the ( "antiferromagnetically" ) 
ordered phase with the exponent (3 = 1/8, the staggered susceptibility diverges with the 
exponent 7 = 7/4. 

Using the elaborated lattice-gas picture we will calculate in the present paper the uniform 
magnetization M(T,h,N) and the uniform susceptibility x{T,h,N) = dM(T,h,N)/dh, 
which are both important and easily accessible quantities for experimental studies. 
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III. UNIFORM MAGNETIZATION AND SUSCEPTIBILITY 



We use Eq. (I2.3P to calculate M(T, h, N) and x(T, h, N) for h > at low temperatures 
T for the frustrated quantum Heisenberg antiferromagnet (12. ip on the considered lattices. 
Using standard formulas for M and \ we get 

M(T, h,N) = -^ [-T In Z(T, h,N)]=M-n, n = T d ^^M) 

where (...) = ^ m=0 x • • • Xlrw=o i(- • -) e ~ /— lgl^ A 4 ' -^0 denotes the grand-canonical 

average for the classical lattice-gas model and n = Ylm=i n ™,- 



A. Frustrated two-leg ladder 

In the one-dimensional case we use the transfer-matrix approach [jjj] to obtain explicitly 
Si g (T, fi, J\f) and hence n and dn/dfi, see Eq. ( 13. ip . After somewhat lengthy but straightfor- 
ward calculations we obtain the following final results in the thermodynamic limit M — > oo 

M(T,h,N) c+ 
N A+' 
T x (T,h,N) d + c\ 



2 ' 



C+ 



1 + w + - w) 2 + 4z 

2 : 
2^ — w(l — w) 

2 V / (l-^) 2 "T4i + 2" ; 



d _ 2z-w(l-2w) _ [2^- W (l-uQ] 2 | w 
+ 2^(1 -w) 2 + 4z 2[(1 -w) 2 + 4z]f 2 

with ^ = e {hi - h)/T , h 1 = J 2 + 2Ji, w = 2e~ Jl/T . 

Some results for M vs h and x vs /i at low temperatures based on Eq. (13.21) are collected 
in Fig. [2] and discussed in Sec. IIII CI 



B. Frustrated bilayer 



The two-dimensional case is more complicated since we do not have an analytical solution 
for the square-lattice lattice-gas model with nearest-neighbor repulsion. Therefore we use 
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FIG. 2: Transfer-matrix approach results for the frustrated two- leg ladder (I2.ip with J\ = 1, 
J% = 4, A/" — > oo: MjM vs /i (upper panel) and Tx/N vs /i (lower panel) at different temperatures 
T = 0.01, 0.1, 0.3. Small filled triangles correspond to the hard-dimer-model predictions valid 
around hi at very low temperatures only. 

the formulas expressing M(T, h, N) and x(T, h, N) in terms of the grand-canonical averages 
n and n 2 , see Eq. (13. ip . and compute n and n 2 using classical Monte Carlo simulations for 
large systems. Namely, we consider systems of M = 50 x 50 and 100 x 100 sites, exploit the 
usual Metropolis algorithm, and perform 220 000 Monte Carlo steps. 

Some results for M vs h and x vs /i at low temperatures based on the Monte Carlo 
simulations of n and n 2 for the lattice-gas model are collected in Fig. [3j We also present 
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FIG. 3: Classical Monte Carlo simulation results for the frustrated bilayer (|2.1h with J\ = 1, 
J 2 = 5, M = 50 x 50 and 100 x 100: MjM vs ft (upper panel) and Tx/M vs ft (lower panel) 
at different temperatures T = 0.01, 0.3, 0.5. Small filled triangles correspond to the hard-square- 
model predictions valid around h\ at very low temperatures only. 

temperature dependences of M and \ at different fields h as they follow from the Monte 
Carlo simulations for the lattice-gas model, see Fig. HI The reported results are discussed 
in Sec. ImCl 
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FIG. 4: Classical Monte Carlo simulation results for the frustrated bilayer (|2.1h with J\ = 1, 
J 2 = 5, M = 50 x 50 and 100 x 100: M/N vs T (upper panels) and T\jM vs T (lower panels) at 
different fields h = 4.0, . . . , 9.1. 

C. Discussion 

Let us discuss the obtained results for the low-temperature magnetization processes for 
the two considered frustrated quantum Heisenberg antiferromagnets. At zero temperature 
T = the magnetization curve M(h) is very simple, namely, it consists only of plateaus 
and jumps, i.e., M{T = 0,h,N)/Af — 1, 1/2 or if h\ < h, h 2 < h < h\ or h < h 2 , cf. 
Ref. 0]. The resulting susceptibility x{T = 0, h, N) is zero in the plateau states, i.e., for all 
values h y£ h 2 and h ^ h 1 . The jumps in M(T = 0, h, N) at h = h 2 and h = hi lead to 5- 
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singularities in the susceptibility. This behavior of the zero-temperature magnetization curve 
reflects the change of the ground states with varying h, i.e., completely empty two sublattices 
for hi < h, one completely occupied with localized magnons sublattice and one completely 
empty sublattice for h 2 < h < hi or completely occupied both sublattices with localized 
magnons for h < h 2 . At h = hi the energy of the (ground) states with different number 
of (independent) localized magnons is the same resulting in a jump in the magnetization 
curve. At h = h 2 the ground-state energy of the states in the subspaces with S z = J\f/2 — r, 
r = 1, . . . , J\f/2 is also the same that provides another jump in the magnetization curve and 
a 5-singularity in the dependence \{T = 0, h, N) vs h at h = h 2 . 

At nonzero temperatures T > the lattice-gas picture implies different scenario for the 
one- dimensional and two-dimensional cases. For the frustrated two-leg ladder the occupation 
of both sublattices with localized magnons increases as h decreases (leading to a decrease of 
M), however, for finite temperatures the two sublattices are equally occupied with localized 
magnons. A smooth decrease of M with two steep parts around hi and h 2 leads to two 
peaks in the susceptibility x a t l° w temperatures, see Fig. EJ 

The situation for the frustrated bilayer is different due to the existence of a long-range 
ordered phase at T/Ji < T /Ji = l/[21n(-\/2 + 1)] « 0.567296 (see the phase diagram 
in Fig. [lb). For a given temperature T < To by variation of the magnetic field (that 
corresponds to a horizontal line in the phase diagram in Fig. [Tb) the system undergoes two 
phase transitions at hi(T) and h 2 (T), where the line T = const, crosses the critical line 
T c (h). [Clearly, hi(T = 0) = hi and h 2 (T = 0) = h 2 .] For the fields between hi(T) and 
h 2 (T) there is a difference in the occupation of the two sublattices, i.e., we have a kind 
of "antiferromagnetic" long-range order of the localized magnons. The total density of the 
localized magnons related to the uniform magnetization M exhibits a kind of inflection, see 
the upper panel in Fig. [31 which results in a sharp peak in x, see the lower panel in Fig. |3j 
An infinite slope in M{h) at hi(T) and h 2 (T) resulting in the 5-peaks in x{h) occurs at 
T = 0, whereas at finite temperatures T the well-pronounced peaks are finite, i.e., there is 
no divergency of x- An alternative way to cross the critical line is to fix h at h 2 < h < hi and 
to vary the temperature T > 0. Corresponding profiles for M and Tx are given in Fig. HI 
Recall that the zero-temperature value of magnetization M(T = 0, h, N)/J\f is 1 if hi < h, 
1/2 if h 2 < h < hi, and if h < h 2 . The susceptibility x is strongly suppressed in the 
"antiferromagnetically" ordered phase. That indicates that the zero-temperature plateau 
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state between hi and h\ extends to finite temperatures (note that this statement does not 
hold for the plateau state at h < hi and h\ < h). Again the temperature profiles of % 
indicate the critical line T c (h), see the lower panels in Fig. |H 

Let us finally briefly illustrate the role of the interacting localized-magnon states by 
comparison with the results considering only independent localized-magnon states, i.e., for 
the lattice-gas models with infinite nearest-neighbor repulsion (the one-dimensional hard- 



dimer and the two-dimensional hard-square model js, 4]). The corresponding data are shown 
in Figs. [2] and [3] by small filled triangles. The hard-dimer and hard-square model predictions 
are correct at low temperatures if h is around hi when the independent localized-magnon 
states do dominate the partition function of the frustrated quantum spin systems. Clearly, 
including interacting localized-magnon states a much wider region of validity is obtained in 
comparison with the hard-core-object description. 

IV. CONCLUSIONS 

To summarize, we have considered the low-temperature magnetization processes for two 
frustrated quantum Heisenberg antiferromagnets. The possibility to examine these systems 
in great detail is owing to emergent Ising degrees of freedom in the strong-coupling regime. 
As a result, the frustrated quantum systems are similar to the classical lattice-gas models or, 
equivalently, the Ising antiferromagnets in a uniform magnetic field. Classical lattice systems 
are much more easier to study either analytically (one-dimensional case) or using classical 
Monte Carlo simulations (two-dimensional case). The most interesting features of the low- 
temperature magnetization processes are related to the finite-temperature order-disorder 
phase transition which occurs in the two-dimensional case. We have found that the field and 
temperature dependences of the uniform magnetization and susceptibility signalize about 
ordering of the localized magnons on vertical bonds of the frustrated bilayer. The jumps and 
the plateaus found at zero temperature for both systems are typical for strongly frustrated 
quantum magnets The steps are smeared out already at quite low temperatures, where 
the plateau is still visible. A prominent feature at finite temperatures are the spikes in 
x{T,h,N) = dM(T, h, N)/dh at the ends of the plateau observed for the ladder system. 
Such spikes have been found theoretically and experimentally, for instance, for the triangular 



spin-1/2 Heisenberg antiferromagnet 15, 16] and for frustrated magnetic molecules [171 ] . For 
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the bilayer model the sharp peaks in x(T, h, N) are related to the phase transition present 
in this system. As a result, they appear for T < T = Ji/[21n(\/2 + 1)] « 0.567 296 Ji, i.e., 
up to quite large temperatures. 

Although in the present paper our consideration is restricted to two particular frustrated 
lattices it can be obviously extended to other lattices, for example, the frustrated three-leg 
ladder in one dimension |6j,[18| or frustrated bilayers consisting of two triangular, honeycomb, 
kagome etc layers in two dimensions. 

In general, bearing in mind the strong-coupling regime for which our consideration is valid, 
we may speak about a system of weakly interacting dimers organized in various lattices 
in different dimensions in such a way that the singlet states on the dimers are localized. 
Frustrated interactions are necessary to achieve this goal. From exp eriment al point of view 



(for solid-state realizations of some similar models see Refs. 19|, |20|) it would be interesting 



to consider the case when the conditions of localization are slightly violated (see, for example, 



Ref. [21] were magnetization curves for such a frustrated two-leg ladder are discussed at 



zero temperature). 
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